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Abstract. We derive the complete asymptotic expansion in terms of powers of N for the 
geodesic /-energy of TV equally spaced points on a rectifiable simple closed curve T in W, p > 2, 
as TV — > co. For / decreasing and convex, such a point configuration minimizes the /-energy 
yijYfr /(^( x j ; x fe)); where d is the geodesic distance (with respect to Y) between points on T. 
Completely monotonic functions, analytic kernel functions, Laurent series, and weighted kernel 
functions / arc studied. Of particular interest are the geodesic Riesz potential l/d s (s ^ 0) and 
the geodesic logarithmic potential log(l/d). By analytic continuation we deduce the expansion 
for all complex values of s. 



1. Introduction 

Throughout this article, T is a Riemannian circle (that is, a rectifiable simple closed curve in 
M. p , p > 2) with length |T| and associated (Lebesgue) arclength measure a — crp. We denote by 
£(x, y) the length of the arc of T from x to y, where x precedes y on T. Thus £(x, y) +£(y, x) = |T| 
for all x, y G r. The geodesic distance d(x, y) between x and y on T is given by the length of the 
shorter arc connecting x and y, that is 



(1.1) d(x,y):=dr(x,y):= m in{£(x,y),£(y,x)} = H 



*(x,y)-E 



The geodesic distance between two points on T can be at most |T|/2. 

Given a lower semicontinuous function / : [0, |r|/2] — > RU{+oo}, the discrete /-energy problem 
is concerned with properties of N point systems z* N , . . . , N on T (N > 2) that minimize the 
/-energy functional 



N N 

(1-2) G / (x 1 ,...,x JV ):=^/(d(x i ,x fc )):=^^/(d(x J -,x fc )), 

j^k j=i k=i 

over all N point configurations u>n of not necessarily distinct points xi, . . . , xjy on T. The following 
result asserts that equally spaced points (with respect to arclength) on T are minimal /-energy 
point configurations for a large class of functions /. 

Proposition 1.1. Let f : [0, |T|/2] ->1U {+00} be a lower semicontinuous function. 

(A) // / is convex and decreasing, then the geodesic f -energy of N points on T attains a 
global minimum at N equally spaced points on T. If f is strictly convex, then these are the only 
configurations that attain a global minimum. 

(B) If f is concave and decreasing, then the geodesic f -energy of N points on T attains a global 
minimum at antipodal systems ojjv with \N /2\ points at p and |_7V"/2J points at q, where p and q 
are any pair of points on T with geodesic distance |T|/2. If f is strictly concave, then these are 
the only configurations that attain a global minimum. 
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Part (A) of Proposition 11.11 follows from a standard "winding number argument" that can be 
traced back to the work of Fejes Toth [33]. The result in the general form stated here appears 
explicitly in the work of M. Gotz [TBI Proposition 9] who uses a similar notion of "orbits." For 
completeness, we present in Section [4] a brief proof of Part (A). 

Remark. Alexander and Stolarsky |2J studied the discrete and continuous energy problem for 
continuous kernel functions / on compact sets. In particular, they established the optimal- 
ity of vertices of a regular _/V-gon circumscribed by a circle C a of radius a for various non- 
Euclidean metrics p(x, y) (including the geodesic metric) with respect to an energy functional 
E a ,\(xi, . . . ,Xjv):=cr([/)(x :; -,Xfc)] A ), < A < 1, on C a where a is an elementary symmetric function 
on (™) real variables. This result does not extend to the complete class of functions in Proposi- 
tion 11.11 and vice versa. However, both cover the generalized sum of geodesic distances problem. 

In the case of Riesz potentials we set 

f s (x):= - x~ s , s<0, / (x):=log(l/a:), f s {x):= X - s , s > 0. 

Then Proposition ll . If A) asserts that equally spaced points are unique (up to translation along the 
simple closed curve P) optimal geodesic f s - energy points for s > — 1. (For s > this fact is also 
proved in the dissertation of S. Borodachov [6j Lemma V.3.1], see also [7].) Proposition 11.1( B) 
shows that for s < — 1 and N > 3, antipodal configurations are optimal /^-energy points, but 
equally spaced points are not. (We remark that if Euclidean distance is used instead of geodesic 
distance, then the A-th roots of unity on the unit circle cease to be optimal /^-energy points when 
s < -2, cf. [5; and [9].) 

For s = — 1 in the geodesic case, equally spaced points are optimal but so are antipodal and 
other configurations. Fejes Toth [15] showed that a configuration on the unit circle is optimal 
with respect to the sum of geodesic distanced (s = —1) if and only if the system is centrally 
symmetric for an even number of points and it is the union of a centrally symmetric set and a set 
{xi, . . . , X2/c+i} such that each half circle determined by Xj (J = 1, . . . , 2k + 1) contains k of the 
points in its interior for an odd number of points. (This result is reproved in [19].) These criteria 
easily carry over to Riemannian circles. In particular, any system of A equally spaced points on 
P and any antipodal system on T satisfy these criteria. 

Remark. Equally spaced points on the unit circle are also universally optimal in the sense of Cohn 
and Kumar [10] . that is, they minimize the energy functional J2j^k /(l x i~ x fc| 2 ) f° r an y completely 
monotonic potential function /; that is, for a function / satisfying (— l) k f^ k \x) > for all integers 
k > and all x £ [0,2]. 

To determine the leading term in the energy asymptotics it is useful to consider the continuous 
energy problem. Let SDT(r) denote the class of Borel probability measures supported on T. The 
geodesic f -energy of \x S 9Jt(r) and the minimum geodesic j '-energy of T are defined, respectively, 
as 

Z?M:= / / /(d(x,y))d At (x)d At (y), Vf(T):=w£ [l g f [n] : fx G 9Jt(r)} . 

The continuous f -energy problem concerns the existence, uniqueness, and characterization of a 
measure /ir satisfying IT) — 2|[/ir]- If such a measure exists, it is called an equilibrium 
measure on T. 

Proposition 1.2. Let f be a Lebesgue integrable lower semicontinuous function on [0, |T|/2] and 
convex and decreasing on (0, |r|/2]. Then the normalized arclength measure Or is an equilibrium 
measure on T and 

(1.3) Jim G } ^P)/N 2 = Vf(T). 

N—KX) 

If in addition, f is strictly decreasing, then err is unique. 



*The analogue problem for the sum of (Euclidean) distances on the unit circle was also studied by Fejes Toth 1141 
who proved that only (rotated copies) of the N-th roots of unity are optimal. 
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The proofs of the propositions in this introduction are given in Section 3] 

Note that (|1.3p provides the first term in the asymptotic expansion of G/(w^P) for large N, that 
is Gf(ujff) ~ Vj!(T) N 2 as N — > oo. The goal of the present paper is to extend this asymptotic 

expansion to an arbitrary number of terms. The case when limAr^oo G f(ujff)/N 2 — > oo as N — > oo 
is also studied. For a certain class of functions / it turns out that the leading term is of the 
form ao2 <^(s ) |r| — s ° N 1+8 ° for some so > lj where ao = lim a ,_ >0 + x s "f(x) is the coefficient of the 
dominant term in the asymptotic expansion of / near the origin and £(s) is the classical Riemann 
zeta function. However, such a leading term might even not exist. Indeed, if the function / 
has an essential singularity at and is otherwise analytic in a sufficiently large annulus centered 
at zero, then the asymptotics of the geodesic /-energy of equally spaced points on L contains an 
infinite series part with rising positive powers of N determined by the principal part of the Laurent 
expansion of / at 0. Consequently, there is no "highest power of N" , see Examples 12.91 and 12.101 
below. 

An outline of our paper is as follows. In Section[2l the geodesic /-energy of equally spaced points 
on r is investigated. In particular, completely monotonic functions, analytic kernel functions, 
Laurent series, and weighted kernel functions / are considered. Illustrative examples complement 
this study. In Section[3l the geodesic logarithmic energy and the geodesic Riesz s-energy of equally 
spaced points on T are studied. The results are compared with their counterparts when d(-, •) is 
replaced by the Euclidean metric. The proofs of the results are given in Section 2] 

2. The geodesic /-energy of equally spaced points on r 

Definition 2.1. Given a kernel function / : [0, |T|/2] — > CU {+oo}, the discrete geodesic f -energy 
of N equally spaced points z^n, . . . , zjy,iv on T is denoted by 

N-l 

A4(T,/;iV):=2/(d(z i ,j V) z M r)) = N £ /(dfz,- ^z^)). 

j/fe 3=1 

Set N = 2M + k [k = 0, 1). Using the fact that the points are equally spaced, it can be easily 
shown that 

LJV/2J 

(2.1) M(T, /; N) = 2N f( n l r l l N ) " (1 " «) /(l r l / 2 ) N - 

71=1 

An essential observation is that the geodesic /-energy has (when expressed in terms of powers of 
N) different asymptotics for even N and odd N. We remark that for real-valued functions / a 
configuration of equally spaced points is optimal with respect to the geodesic /-energy defined in 
(|1.2p . whenever / satisfies the hypotheses of Proposition 1 1 . 1 { A) . 

An application of the generalized Euler-MacLaurin summation formula (see Proposition 14.11 
below) yields an exact formula for A4(T, /; N) in terms of powers of N. The asymptotic analysis 
of this expression motivates the following definition. 

Definition 2.2. A function / : [0, |T|/2] -^CU {+00} is called admissible if the following holds: 

(i) / has a continuous derivative of order 2p + 1 on the interval (0, |T|/2]; 

(ii) there exists a function S q (x) of the form S q (x) = X)n=o a ™ x~ Sn , where a n and s n (n = 
0, . . . , q) are complex numbers with Re sq > Re s\ > ■ ■ ■ > Re s q Q and Re s q + 2p > or 
s q = — 2p such that for some S > 

(a) 1 - Res g + S > 0, 

(b) f {7(2/) - S q (y)} dy = 0(x 1+s -^ ) as x -> 0+, 
Jo 

(c) {/(x) - S q (x)} {u) = 0{x s - s "- v ) as x -> 0+ for all v = 0, 1, . . . , 2p + 1. 



tThe powers in S q (x) are principal values. 
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For p > 1 an integer the following sum arises in the main theorems describing the asymptotics 
of M(TJ;N): 

(2.2) B p (T, /; N)-~N 2 £ ^jff (|r| /AT) 2 " / (2 "- 1} (|r| /2), JV = 2M + «, « = 0, 1, 
where B m (x) denotes the Bernoulli polynomial of degree m defined by 



-e 



ml 

m=0 k=0 



where Bq = 1, B\ = —1/2, ...are the so-called Bernoulli numbers. Recall that B2k+i = 0, 
(-l)*- 1 ^* > for k= 1,2,3,..., and B n (l/2) = (2 1 "™ - l)B n for n > (pQ). 

Theorem 2.3 (general case). Let f be admissible in the sense of Definition \2.2\ and suppose none 
of so, si, . . . , s q equals 1. Then, for N — 2M + k with k = or k = 1, 

(2.3) M(T, /; N) = V f (T) N 2 + £ a„^^iV 1+s " + B p (T, /; AT) + fl^r, /; N), 

n=0 ' ' 



where 

(2-4) V f (T) = pE a " Ui r^ ST , + ]f| / ' " if-S q )(x)dx 



2 « (in/2) 1 - 5 " 2 Hri/2 



n=0 



and the remainder term satisfies 9ip(T,f;N) = 0{N 1 ~ 2p ) + 0(N 1 ~ S+Sq ) as N -> oo if 2p ^ 
<5-Res g; w/iereas «H p (r,/;iV) = ©(iV 1 "^ log iV) i/2p = <J - Res 9 . 

The next result involves the Euler-Mascheroni constant defined by 

,. (, 1 1 1 1 , 

7:= hm 1 + - + - + - H 1 logn 

n— too \ 2 3 4 n 

Theorem 2.4 (exceptional case). Lei / 6e admissible in the sense of Definition [2~B and s g / = 1 
for some 1 < g' < gll Then, for N — 2AI + k with k = or K = 1, 

M(T, /; TV) = A v at2 log at + V)(T) iV 2 + ^ a „^^Ar 1+s » + B p (r, /; N) + m p (T, /; N), 

' ' n=0, ' ' 

where 

(2.5) ^f(r) = p| E a » [ _ Sn + X (/-S 9 )(z)&-<v(log2- 7 ) j 

and i/ie remainder term satisfies V\ p (T,f;N) = OiN 1 -^) + 0(N 1 ~ S+S «) as N -> 00 if 2p ^ 
5-Rcs ? , whereas 9\ P (T, /; AT) = <D(N 1 ~ 2p log A^) i/ 2p = <5 - Re s 8 . 



Remark. Both Theorems 12 .31 and 12 .41 show that only the coefficients of the nonpositive even powers 
of AT depend on the parity of N. These dependencies appear in the sum B P (F, /; N). 

Remark. If f(z) = S q (z) = X)n=o an2:_S " ^ or some Q an d Reso > ■ • • > Res g , then all expressions 
in Theorems 12.31 and 12.41 containing f — S q vanish. In general, the remainder term 9t p (r, /; N) is 
of order 0(N 1 ~ 2p ), where the integer p satisfies Res g + 2p > 0. In particular, this holds for the 
Riesz kernels (cf. Theorems 13.31 and 13.51 below) . 



'By Definition 12. 21 there is only one such s q i 
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Completely monotonic functions. A non-constant completely monotonic function / : (0, oo) — > 
R has derivatives of all orders and satisfies (— l) k f^(x) > (cf. [12])@ In particular, it is a con- 
tinuous strictly decreasing convex function. Therefore, by Proposition ll.il equally spaced points 
are optimal /-energy configurations on the Riemannian circle T. 

By Bernstein's theorem [29, p. 161] a function is completely monotonic on (0,oo) if and only if 
it is the Laplace transformation f(x) — J Q e~ xt dfi(t) of some nonnegative measure /i on [0,oo) 
such that the integral converges for all x > 0. 

The following result applies in particular to completely monotonic functions. 

Theorem 2.5. Let f be the Laplace transform f(x) — e~ xt d/i(t) for some signed Borel mea- 
sure /i on [0, oo) such that J °° t m d\/j,\(t), m — 0, 1, 2, ... , are all finite. Then for all integers p > 1 
and N = 2M + k with n = 0, 1 

M(T, /; AT) = {A £ 1 ~ e ^' r ' /2 Mt)} N^(-ir^ 2 -^N^+B p (T, /; N)+0(N^), 
where n m '-= J t m dfj,(t) denotes the m-th moment of fi. 

Remark. The derivation of the (complete) asymptotic expansion for A4(T, f; N) as N — > oo for 
Laplace transforms for which not all moments fj, m are finite, depends on more detailed knowledge of 
the behavior of f(x) near the origin. For example, for integral transforms G(x) — J h(xt)g(t) dt 
there is a well-established theory of the asymptotic expansion of G(x) at + . See, [T7], [H], [4] or 
[24] and [13]. These expansions give rise to results similar to our theorem above. 

Remark. Recently, Koumandos and Pedersen [21] studied so-called completely monotonic functions 
of integer order r > 0, that is functions / for which x r f{x) is completely monotonic. The com- 
pletely monotonic functions of order are the classical completely monotonic functions; those of 
order 1 are the so-called strongly completely monotonic functions satisfying that (— l) k x k+1 f^ k \x) 
is nonnegative and decreasing on (0,oo). In [21] it is shown that / is completely monotonic of 
order a > (a real) if and only if / is the Laplace transformation of a fractional integral of a 
positive Radon measure on [0, oo); that is 

f°° 1 /"* 

f(x)= e~ xt J a {n}(t)dt, J a [ IJt \(t):=—- {t- S ) a - 1 dix{ S ). 



T(a) 

Results similar to Theorem 12 . 5 1 hold for these kinds of functions. However, the problem of giving 
an asymptotic expansion of f(x) near the origin is more subtle. 

Analytic kernel functions. If / is analytic in a disc with radius |r|/2 + e (e > 0) centered at 
the origin, then / is admissible in the sense of Definition 1 2 . 21 and we have the following result. 



Theorem 2.6. Let f(z) = X^o a ™ z ™ be analytic in \z\ < |r|/2 + e, e > 0. Then for N = 2M+K 
with k = or k = 1 

M(T, f; N) = | p jf f{x) dx\N* + £ a „ J^AT 1 -" + B P (T, /; N) + O p . ]TU {N^)- 

Note that C(0) = -1/2 and C(~2fc) = for k = 1,2,3, ... . 
Example 2.7. If f(x) — e~ x , then for any positive integer p: 

M(T f ■ N) = — ( 1 - e-l r l/ 2N l N 2 - N + V - 2 C(l-2n) 2 _ 2n 

(2n)! { r\^ N + °»nA N ) 

as N = 2M + k — > oo, where the notation of the last term indicates that the O-constant depends 
on p, \T\ and /. Since f(x) is a strictly decreasing convex function, by Proposition 1 1 . If A) . equally 



§ A completely monotonic function on (0, oo) is necessarily analytic in the positive half-plane ( 29 ). 
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spaced points are also optimal /-energy points. Thus, the relation above gives the complete 
asymptotics for the optimal iV-point geodesic e - (')-energy on Riemannian circles. 

Laurent series kernels. If f(z) is analytic in the annulus < \z\ < \T\/2 + e (e > 0) with a pole 
at z = 0, then / is admissible in the sense of Definition 12.21 and we obtain the following result. 



Theorem 2.8. Let f be analytic in the annulus < \z\ < \T\/2 + e (e > 0) having there the 
Laurent series expansion f(z) — XmL-ic a « z ™> K >\. 

(i) // the residue a_i = 0, then for N = 2M + k with k = 0, 1 

2 P 7(~(-n) 

M(TJ;N) = V f (T)N 2 + ]T a n J±—l N ^ + Bp (T, f; N) + O p ,\ r \ jiN 1 ^), 

n=-K, I I 

where the N 2 -coefficient is 

2 ^ (|r|/2) 1+n 



n 

■ ■ n=-K 

(ii) If the residue a_i ^ 0, then for N = 2M + k with k = 0, 1 

irr" 

where the N 2 -coefficient is 



2p 

M(T, /; IV) = ,!„_, «» lo gW + y, (r) IV 2 + £ a „Bifcj> A"-' + B,(r, /; IV) + O r , n ,i.N 1 -*) 



n#-l' 

Next, we give two examples of kernels / each having an essential singularity at 0. Such kernels 
can also be treated in the given framework, since they satisfy an extended version of Definition 12. 2\ 
see Proof of Examples 12.91 and 12.101 in Section |U 

Example 2.9. Let f(x) = e}l x = Y^Lo l/(nlx n ), x E (0, +oo), /(0) = +oo. We define the entire 
function 

F(z):= V ^L z n = - 7 z- — I e z/w ^(l-w)dw, zeC, 
^ nl 2™J\ W \= P <1 

where ip(z) denotes the digamma function and we observe that, because of < ((n) — 1 < c2~" 
for all integers n > 2 for some c > 0, 

Fix) = e x - 1 - x + V 1 x n = e x + 0(e x/2 ) as x -> oo. 

, n! 

n— 2 

Then 

M(T, /; iV) = 2NF{N/ \T\) + -^N 2 \ogN + V f {T) N 2 - N 

+ £ ( ^ffL ^'- 2m / (2m - 1) (|r| /2) + Ptl r U (N^), 

where 

1 1 n=2 1 1 



e 2/|r| {l-2 7 + log |r|+Ei(2/|r|)}, 
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where Ei(ir) = — e H 1 dt is the exponential integral (taking the Cauchy principal value of 
the integral). In particular it follows that 

ton M{T ±V = 2. 

Since / is a strictly decreasing convex function on (0, oo), by Proposition 1 1 . 1 \ A) . equally spaced 
points are also optimal. Thus, the above expansion gives the asymptotics of the optimal TV-point 
energy. 



e V(0 



Example 2.10. Let Jfc(A) = (— l) fe J_fc(A):=^ cos(k8 ~ \ sin 9) d9 denote the Bess el function 
of the first kind of order k whose generating function relation is given by (cf. [27[ Exercise 5.5(10)]) 



f(x) = exp 



- \ x 



= J2 J n(A>™ for \x\ > 0. 



71— — OO 



For integers m > 2 we define the entire functions 

oo oo oo 

F m {z):= J-n(A) ({n)z n = E G m (z/k), G m {z):= E J-»(A)* n , z e C. 



k=l 



If A is a zero of the Bessel function J_i, then for positive integers p and m > 2 there holds 

m— 1 

M(r, /; zv) - 2iv^ m (/v/ |r|) + 2 E J-n(A) CW |rp" tv 1+ " + v> (r) tv 2 + |r| B 2 (?)/'(|r| /2) 



E 



where 



25 2 „/ 2 "- 1 (|r|/2) 

2n (2n-l)! 
Vjf(r) 



+ 2 J 2 „_i(A) C(l - 2n) ^ irl 2 "" 1 /V 2 - 2 " + ©(/V 1 -^) 



^— ' 1 + n 



n— — oo, 
n#±l 



If, in addition, A < 0, then f(x) is a strictly decreasing convex function and, therefore, M(T, /; TV) 
is also the minimal TV-point /-energy on T and it follows from the observation 



G m [x/k) = exp 





fx 


k 


4 


\k 


X 



that 



lim 



- E MX){-x/k) n , k = 1,2,3,. 

n— — oo 

M(T,f;N) n 



w^co TV/(-TV/|r|) 

If A is not a zero of J_i, then the above asymptotics must be modified to include a logarithmic 
term. 

The weighted kernel function ff{x) — x~ s w(x). Given a weight function w(x), the kernel 
fT( x ) = x~ s w(x) gives rise to the so-called geodesic weighted Riesz s-energy of an TV-point con- 
figuration (xi, . . . , xjv) 

w(d(xj,x k )) 



G™(x 1 , 



,xjv 



[d(xj-,x fc )]' 



For the Euclidean metric the related weighted energy functionals are studied in [5]. 

If w{x) is such that ff(x) is admissible in the sense of Definition 12. 2i then Theorems 12 . 31 and 12.41 
provide asymptotic expansions for the weighted geodesic Riesz s-energy of equally spaced points 
on a Riemannian circle T, which are also optimal configurations if ff{x) is strictly decreasing and 
convex (cf. Proposition II. 1( A)). 
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Theorem 2.11. Letw(z) = E,T=o a « z " be analytic in \z\ < |r|/2+e, e > 0. Set ff(z):=z~ s w{z). 
Then for integers p,q > and s £ C, s not an integer, such that q — 2p < Re s < 2 + q we have 

M(T, ff- N) = V fr (r) N 2 + J2 an 2 ^~ n n V +'-" + B p (T, f?;N) + ^(T, f?;N), 

where B p is defined in (|2.2[) . TTie N 2 -coefficient is the meromorphic continuation to C o/ i/ie 

ri /2 

geodesic ff -energy ofT given by (2/|r|) ff{x)dx for < s < 1; £/iai is 

2 A (|r|/2) 1+ — 

%( r )-p2. a " 1 + n _ s » ^1,2,3,.... 

1 1 n— 

T/ie remainder m p (T, f™;N) is of order 0{N 1 -^) + 0(N s - 2 p) as N -> oo. 

Remark. For s is a positive integer the series X^^Lo a nZ n ~ s is the Laurent expansion of /(z) in 
< |z| < |r|/2+e and Theorem l2. 81 applies. For s is a non-positive integer the series J2^=o a nZ n ~ s 
is the power series expansion of f(z) in < \z\ < \T\/2 + e and Theorem 12.61 applies . 

Example 2.12. Let w(z) — sin(z7r/|L|). Then for Res > not an integer 

°° f_-\\n 

fT(z) = x~'w{z) = £ T^ry (-/ |r|) 2n+1 z 2 " +1 - s 

and, by Theorem 12.111 the geodesic weighted Riesz s-energy of N equally spaced points has the 
asymptotic expansion (0 < Re s < 1 + 2p) 

M(T, f?;N) = Vfw (L) N 2 +(n/ \T\) S £ ff- l)! ^ ^it*-^ N 2+S ~ 2k +B p (r, ff- N)+m p (T, f?;N), 

where B P (T, ff;N) is given in (|23|) , The remainder m p {T, ff; N) is of order 0(iV 1 - 2 P)+e>(iV s - 2p ) 
as N — > oo and 

v m ^ /^f- (-i) fc - x (V2) 2fe - 

For < s < 1 we have 

VfT (r) - |pj ^ /. (*) - - 2 s lF2 ^ _ s/2; 3/JJI - (tt/4) 

expressed in terms of a generalized iF2(; )-hypergeometric function, which is analytic at s not an 
even integer. Hence, Vfw (T) is the meromorphic continuation to the complex plane of the integral 

ITT /(J r|/2 ' ffi x ) dx - We observe that for s = 1/2 we have V fs >{T) = 2,/2/\r\S(l), where S(u) is 
the Fresnel integral S(ii):= JJ* sin(x 2 7r/2) dx. 

As an application of the theorems of this section, we recover results recently given in [5] regard- 
ing the complete asymptotic expansion of the Euclidean Riesz s-energy C s (N) of the JV-th roots of 
unity on the unit circle S 1 in the complex plane C. Indeed, if \z— w\ denotes the Euclidean distance 
between two points £ and z in C, then from the identities \z — £| 2 = 2(1 — cos^) = 4[sin(?/>/2)] 2 , 
where ip denotes the angle "between" £ and z on S , we obtain the following relation between 
Euclidean and geodesic Riesz s-kernel: 

d(C,*) 



|z-Cr s = |2(l-cos^)| s/2 = 
Thus, for (,z € S 1 there holds 



ib 

2 sin- 



2 sin ■ 



C.2 6S 1 . 



(2.6) |z-Cr S = /r(d(C,^)), ^(x):=(sinc|) ff (x) = x~ s smc~ s (x/2), 

where the "sine" function, defined as sincz = (sinz)/z is an entire function that is non-zero for 
\z\ < 7r and hence, has a logarithm g(z) = log sine z that is analytic for \z\ < n (we choose the 
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branch such that logsincO = 0). The function sine s (z/2):= exp[— slogsinc(z/2)] is even and 
analytic on the unit disc \z\ < 2ir and thus has a power series representation of the form 

oo 

sin<T a (z/2) = a n(s)z 2n , \z\ < 2tt, s E C. 

n—0 

It can be easily seen that for s > —1 and s / the function (sgns)f™(x) is a convex and 
decreasing function. Hence, application of Proposition II. If A) reproves the well-known fact that 
the N-th roots of unity and their rotated copies are the only optimal /^-energy configurations 
for s in the range (—1, 0) U (0, oo). (We remind the reader that, in contrast to the geodesic case, 
in the Euclidean case the AT-th roots of unity are optimal for s > —2, s ^ 0, and they are 
unique up to rotation for s > — 2, see discussion in [§].) The complete asymptotic expansion of 
C S (N) — A^S 1 , /™; N) can be obtained from Theorem l2.11l if s is not an integer, from Theorem l2.8l 
if s is a positive integer, and from Theorem 12.61 if s is a negative integer. (We leave the details 
to the reader.) For s £ C with s ^ 0,1,3,5,... and q — 2p < Re s < 2 + q, the Euclidean Riesz 
s-energy for the N-th roots of unity is given by (cf. [H Theorem 1.1]) 

v ' n— 1 v ' 



(2.7) C S (N) = V S N 2 + ^P^N 1+S + Y a n (s) ;}\ ' iV 1+s ' 2 " + ©(iV 1 "^) + 0(N s ~ 2 p) 



as N —> oo, where (cf. [9]) 

r98 , v 2-r((l-.s)/2) (-i)"^(V2) 

Here, B^\x) denotes the generalized Bernoulli polynomial, where B n {x) — S„ (x). Notice the 
absence of the term B P (T, ff] N), which follows from the fact that odd derivatives of ff{x) evalu- 
ated at 7r assume the value 0. (This can be seen, for example, from Faa di Bruno's differentiation 
formula.) 

The entirety of positive odd integers s constitutes the class of exceptional cases regarding 
the Euclidean Riesz s-energy of the iV-th roots of unity. For such s Theorem I2.8( ii) provides 
the asymptotic expansion of C S {N) = ^(S 1 , ff; N), which features an N 2 log TV term as leading 
term. That is, for s = 2L+1, L = 0, 1, 2, . . . , we have from Theorem E^ii) that (cf. !9, Thm. 1.2]) 

(2.9) C S (N) = ^M N 2 logN + V fr (S 1 )N 2 +y2 a m ( s ) 2 -^-^N 1+s ~ 2m + OiN 1 ^), 

11 m=Q, O) 

where the coefficients a m (s) are given in (|2.8p and 

%( §1 ) = ^ E a ^ s U'L , i -M S )(log2- 7 ) 




We remark that in [9, Thm. 1.2] we also give a computationally more accessible representation of 
Vf-w (S 1 ). The appearance of the N 2 log N terms can be understood on observing that the constant 
V s in (|2.7p has its simple poles at positive odd integers s and when using a limit process as s — > K 
(K a positive odd integer) in l|2.7[) . the simple pole at s = K need to be compensated by the 
simple pole of the Riemann zeta function in the coefficient of an appropriate lower-order term. 
This interplay produces eventually the N 2 log N term. 

3. The geodesic Riesz s-energy of equally spaced points on r 

Here, we state theorems concerning the geodesic Riesz s-energy of equally spaced points on 
r that follow of the results from the preceding section together with asymptotic properties of 
generalized harmonic numbers. The proofs are given in Section 2J 



'The function sgns denotes the sign of s. It is defined to be —1 if s < 0, if s = 0, and 1 if s > 0. 
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Definition 3.1. The discrete geodesic Riesz s- energy of N equally spaced points z^jv, . . . , z^at 
on r is given by 

N-l 

M S {T; N):= £ [d(z jfN , z kiN )}~ s = N ^ n, *n,n)]~' , a G C. 

The discrete logarithmic geodesic energy of N equally spaced points z^at, . . . , zjv,N on r enters in 
a natural way by taking the limit 

M S (T;N)-N(N-1) 1 



(3.1) Mo g (r^):=lim^i^^l_iZ =Elog 

s f-f d(z J -.7v,Zfc ] Ar) 

We are interested in the asymptotics of A4 S (T; N) for large iV for all values of s in the complex 
plane and we shall compare them with the related asymptotics for the Euclidean case given in our 
recent paper [9] . In the following we use the notation 

z»= / / v s °(ry.=M{xiM : /x g 9Jt(r)}, 

3.1. The geodesic logarithmic energy. 

Theorem 3.2. Let q be a positive integer. For N = 2M + K, K = 0, 1 



M og (r ; AO = v* g (r) iv 2 - iviogiv + iviog g - £ (ffz^ 22 "^ 2 " 2 " + <V(^~ 2<? ) 

as iV -> oo. ffere, V£ g (T) = 1 - log(|r|/2). 

Remark. The parity of AT affects the coefficients of the powers jV 2_2m , m > 1. The A^ 2 -term 
vanishes for curves T with |T| = 2e and the TV-term vanishes when |T| = 2tt. By contrast, the 
Euclidean logarithmic energy of N equally spaced points on the unit circle is given by (cf. [9]) 

C log (N) = -N log N. 

3.2. The geodesic Riesz s-energy. The next result provides the complete asymptotic formula 
for all s ^ 1. This exceptional case, in which a logarithmic term arises, is described in Theorem l3.5l 

Theorem 3.3 (general case). Let q be a positive integer. Then for all s G C with s / 1 and 

Re s + 2q > there holds 

(3.2) 

M S (T; N) = V'F) N 2 + ^N 1 * - £ {s) 2n _^N^ + 0^^) 

as N -> oo, where Kf(r) = (|r|/2)~ s /(l - s) and N = 2M + k, k = 0, 1. 

In (13.21) the symbol (s) n denotes the Pochhammer symbol defined as (s) = 1 and (s) n+1 = 
(n + s)(s) for integers n > 0. 

Remark. It is interesting to compare (|3.2[) with (|2.7p . It should be noted that in both the geodesic 
and the Euclidean case, the respective asymptotics have an A^-term whose coefficient is the 
respective energy integral of the limit distribution (which is the normalized arc-length measure) or 
its appropriate analytic continuation, and an A 7 " 1+s -term with the coefficient 2 £(s)/|r| s . Regarding 
the latter, it has been shown in [55] that for s > 1 the dominant term of the asymptotics for the 
(Euclidean) Riesz s-energy of optimal energy Appoint systems for any one-dimensional rectifiable 
curves in R p is given by 2 £(s)/|r| s iV 1+s . Regarding the remaining terms of the asymptotics of 
/A S (T;N) and C S (N) one sees that the exponents of the powers of N do not depend on s in the 
geodesic case but do depend on s in the Euclidean case. 
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Remark. In the general case s 7^ 1, the asymptotic series expansion (|3.2p is not convergent, except 
for s = 0, —1, —2, . . . when the infinite series reduces to a finite sum. The former follows from 
properties of the Bernoulli numbers and the latter from properties of the Pochhammer symbol 

For a negative integer s we have the following result. 
Proposition 3.4. Let p be a positive integer. Then 



2n nj2 — 2n 



p + 1 p + 1 ^ V 2n j 

n— 1 v 7 

+ ffl( J B p+1 ( K /2)-B p+1 )7Vi-P 

J5 + 1 

/or N = 2M + k, k = 0, 1. T/ie right-most term above vanishes for even p. 

Remark. The corresponding Euclidean Riesz (— p)-energy of iV-th roots of unity reduces to 

£- p (N) = V- P N 2 if p= 2,4,6,.... 

Remark. The quantity Al_i(§; TV) gives the maximum sum of geodesic distances on the unit circle. 
Corollary 13.41 yields 



(3.3) M-i{S;N) = ^(N 2 - k) , N = 2M + k, k = 0, 1. 

We remark that L. Fejes Toth [15) conjectured (and proved for N < 6) that the maximum sum of 
geodesic distances on the unit sphere S 2 in R 3 is also given by the right-hand side in (|3.3j) . This 
conjecture was proved by Sperling [55] for even iV0 and by Larcher g3] for odd iV0 An essential 
observation is that the sum of geodesic distances does not change if a given pair of antipodal 
points (x, x') is rotated simultaneously, since d(x, y) + d(x',y) = n for every y G § 2 . 

In the exceptional case s = 1 a logarithmic term appears. 

Theorem 3.5. Let q > 1 be an integer. For N = 2M + k, k = 0, 1, 

Ml (T; N) = log TV — l0g '7 7 7V 2 -If B2 " (K/2) 2 2 "iV 2 - 2 " 

(3 .4) > r l l r l/ 2 l r l^ 2n 

2 B 2q +2(K/2) 2 +2 2 

-W|f| 2g + 2 2 * , 

where < q ,Ar, K < 1 depends on q, N and k. 

Remark. A comparison of the asymptotics (|3.4p and the corresponding result for the Euclidean 
Riesz 1-energy of iV-th roots of unity (cf. ([13]) and [9, Thm. 1.2]), 

A(A0 = iiV 2 log N + 7 " lQS(7r/2 V 2 + V (-l)"^(l/2)2C(l-2n) jv2 _ an + Q 
7T 7T ^ (2n)! (2tt) 

shows that for |T| = 2ir the dominant term is the same and the coefficients of all other powers of N 
differ. The latter is obvious for the iV 2 -term, and for the 7V 2_2 ™-term, follows from the fact that the 
coefficient in (]3.4p multiplied by it is rational whereas the coefficient in the asymptotics for C%(N) 
multiplied by 7r is transcendental. Interestingly, except for s = 1, there are no other exceptional 
cases with an N 2 logiV term in the asymptotics of A4 S (T; N), whereas in the asymptotics of C S (N) 
there appears an N 2 logiV term whenever s is a positive integer, cf. [HI Thm. 1.2]. 



"Sperling mentions that his proof can be easily generalized to higher-dimensional spheres. 
"Larcher also characterizes all optimal configurations. 
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4. Proofs 

Proof of Proposition ] Part (A). The proof utilizes the "winding number" argument of L. 
Fejes Toth. The key idea is to regroup the terms in the sum in (|1.2p with respect to its to nearest 
neighbors (to = 1, . . . , N) and then use convexity and Jensen's inequality. 

W.l.o.g. we assume that wi, . . . ,Wjy on T are ordered such that w*, precedes Wk+i (denoted 
Wfc -< Wfc + i). We identify mvj+n with Wj for j = 1, . . . , N — 1. By convexity 

JV JV JV-1 N N-l N 

j=l k=l k=l j=l k=l j=l 

Let zi t jv -<•••-< zjv.at be N equally spaced (with respect to the metric d) points on T. Set 
zq,jv = zn.n- Assuming further that this metric d also satisfies 



(4.2) i^d(x i ,Xj +fc ) < d(z ,Ar,z fe ^), k = 1, . . . , N - 1, 

i=i 

for every ordered iV-point configuration xi -< •• • -< xjv with x 3 = x^+jv, it follows that 

N-l 

G / (w 1 ,...,w JV )>iV^/(d(z 

fc=i 

It remains to show that the geodesic distance satisfies (|4.2p . From 

d(xj, x fc ) = min{^(x 3 , x fc ), |r| - £(Xj, x fc )} if < fc - j < N 
and additivity of the distance function £(■, •) it follows that 



f N N k 

N 



^d( Xj ,X j+fc ) < 



E ^( X i' X J+fe) = E E ^( x J+n-l> x J+n) = l r l fc ' 
j=l j=l n=l 

^(|r|-^(x i) x J - +fe )) = |r|(iv_fc) 



and therefore 

1 - 

-^d( Xj ,x J+fc ) < min{|r|fc/7V,|r| (TV - fc)/iV} = d(z ,A r ,z fe , J v). 
j=l 

In the case of a strictly convex function / we have equality in (|4.ip if and only if the points 
are equally spaced. This shows uniqueness (up to translation along the simple closed curve T) of 
equally spaced points. 

Part (B). Given N = 2M + k (k = 0, 1) let wjv denote the antipodal set with M + k points 
placed at the North Pole and M points at the South Pole of T, where both Poles can be any two 
points on T with geodesic distance |T|/2. Thus, the geodesic distance between two points in lon 
is either or \T\/2. Hence 

(4.3) G f (io N ) = 2M (M + k) f(\T\ /2) = i/(|r| /2) (N 2 - «) . 

Since adding a constant to G/ does not change the positions of optimal /-energy points, we 
may assume w.l.o.g. that /(0) = 0. In fact, we will prove the equivalent assertion that if / 
is a non-constant convex and increasing function with f(0) = 0, then the functional Gf has a 
maximum at con, which is unique (up to translation along T) if / is strictly increasing. (Note that 
by these assumptions f(x) > 0.) Indeed, any iV-point system Xm of points Xi, . . . ,xat from P 
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satisfies 

Gf(x N ) = / ( in 12) e f{ f^ ] < m\ /2) £ ^f 1 - m w q ^jr 

</(|r|/2)^|^ = i/(|r|/2) (n 2 - k ), 

where we used that antipodal configurations are optimal for the "sum of distance function" (/ is 
the identity function id) and relation (|4.3p with / = id. Note that the first inequality is strict if 
there is at least one pair (j, k) such that < d(xj,Xfc) < \T\/2. On the other hand, if Xn — lon, 
then equality holds everywhere. □ 

Proof of Proposition For Lebesgue integrable functions / the minimum geodesic /-energy 
V f 9 (T) is finite, since 2* [or] = J/(d(x,y)) dor(x) = (2/|r|) jf l/2 f(£) d£ ^ oo (y € T arbitrary). 
Moreover, for lower semicontinuous functions /, a standard argument (see [22) ) shows that the 
sequence {Gf{uj^ )/[N(N — 1)]}at>2 is monotonically increasing. Since / is Lebesgue integrable, 
this sequence is bounded from above by Jj[ar]\ thus, the limit limjv^oo Gf(Loff)/N 2 exists in 

this case. If / also satisfies the hypotheses of Proposition 1 1 . If A) . then liniAr^oo Gf(uj^)/N 2 = 
Tj[<Jr]- (By a standard argument, one constructs a family of continuous functions F £ (x) with 
F £ (x) = f(x) outside of e-neighborhoods at points of discontinuity of /, f{x) > F £ (x) everywhere 
and lim 6 _+o F £ (x) = f[x) wherever / is continuous at x. Then the lower bound follows from 
weak-star convergence of ^[w^P] as N — > oo and, subsequently, letting e — > 0.) □ 

We next present some auxiliary results that are needed to prove the main Theorems 12 .31 and 12 .41 
We begin with the following generalized Euler-MacLaurin summation formula. 

Proposition 4.1. Let lo = or lo — 1/2. Let M > 2. Then for any function h with continuous 
derivative of order 2p + 1 on the interval [1 — lo, M + lo] we have 

£>(*) = I h(x)dx + (1/2 -c,){h(a) + h(b)} + ^ f^{h (2k - 1 Hb)-h^- 1 Ha)} 



k=l ua fe=l 

b 

C 2p+l {x)h {2p+l \x)dx, a=l-w,b = M- 



1 



(2p + l)\ 

where Ck{x) is the periodized Bernoulli polynomial Bk(x — [x\). 

Proof. For lo = 0, the above formula is the classical Euler-MacLaurin summation formula (cf., 
for example, [3]). For lo = 1/2, iterated application of integration by parts yields the desired 
result. □ 

Let / have a continuous derivative of order 2p + 1 on the interval (0, |T|/2]. Then applying 
Proposition O with h(x) = f(x\T\/N) and w = k/2, where N = 2M + k > 2, K = 0, 1, we obtain 

LiV/2j N/2 

M(T, /; N) = 2N £ /(« \?\ /*) "(!"«) /(|r| /2)N = 2N / f(x\T\/N) dx 
n=l Jl - U 



- - J] N{f((l - lo)\T\/N) + /(|r|/2)} + 2iV £ {f(x\T\/N)Y 2k - V 

k=l 

2N f N/2 „ ,,,.,.„ r 2 „+ii , , . _ fl 



1 I tr ^ 

C 2p+1 (x) {f(x\T\/N)} {2p+1) [x)dx-2( l -~Lo\ /(|r| /2)N. 



(2p+l) 



Regrouping the terms in the last relation and using the fact that B 2 k+i = B 2 k+i(l/2) = for 
k = 1, 2, 3, . . . and B\{lo) = lo — 1/2, we derive the exact representation 



2 



|r|/2 



(4.4) M(T, f;N)=N Z j- f(y) dy - A P (T, /; N) + B P {T, /; N) + 72, (T, /; N) 

I 1 1 J(i-u)\r\/N 
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valid for every integer N > 2, where 

A P (TJ;N):= - 2B 1 (u>)f((l-uj)\T\/N) - 2Nj2^nT {/Nrl/iV)}^ 



2p 



(4.5a) 



= i\N 2 £ ^P- (iri W - -) |r| /N), 

E^flri/^/^dri/^ 



(4.5b) B p (TJ;N):~ N ^ ^ 



fc= 



(4.5c) ft p (r,/;iV):=27V 



(|r|/iv) 



2p+l .JV/2 



^p+i^/^+^xiri/AO^. 



(2p + l)! A- 

If / is admissible in the sense of Definition 12.21 then by linearity 

M(T, f;N)=M (r, S q ; N) + M (T, f - S q ; N), 

where the term A4(T, S q ;N) contains the asymptotic expansion of M.(T, f; N) and the term 
A4(T,f — S q ;N) is part of the remainder term. The next lemma provides estimates for the 
contributions to the remainder term in the asymptotic expansion of M(T, /; N) as iV — > oo. 

Lemma 4.2. Let f be admissible in the sense of Definition 1 2. 2\ Then as N — > oo: 



(i-uO|r|/jv 



(f-S q )(y)dy = 0(N l - s+s «), 

A p (TJ-S q ;N)=0(N 1 - s +^), 
'0{N l ~ 2p ) 



n p {TJ-S q] N) 



The O-term depends on \T\, p, s q , and f. 



if2p^S — Res q , 
OiN 1 - 2 * log N) if2p = 5-Res q . 



Proof. The first relation follows directly from Definition l2.2f ii.aV The second estimate follows 
from Definition I2.2f ii.b) and (|4.5a[) ; that is for some positive constant C 

2p 

\A P (T, f-S q ;N)\< ^-N 2 £ (|r| /N) r (/ - S q )^((l u ) \T\ /N) 



T — n 

1 1 r=l 



2p 



\B r (w) 



(|r| /N) r (1 - (|r| /N y+s-R C s q - r +i 



The last estimate follows from Definition 12. 2f ii.b). (|4.5c[) and the fact that 
(4.6) \C 2p +i{x)\ < (2p+ 1) \B 2p \ for all real x and all p = 1, 2, . . . ; 

that is for some positive constant C 



\K p {Y,f-S Q ;N)\ <2N 



(\T\/N) 



2p+l f N/2 



(2p + l)\ a_ & 

B 2p 
(2p)\ 



<2CN^-{\T\/N) S - Rcs 



C 2p+1 (x)\ (f-S q )^ +1 \x\T\/N) 

N/2 



dx 



1— u 



□ 



Other functions arising in the asymptotics of A4(T, /; N) are defined next. 
Definition 4.3. Let oj = 0, 1/2 and p be a positive integer. For s 6 C with s / 1 



2p 



, 1 - ;• 



r=0 



(a) 2P +i r 

(2p+l)l 
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which we call incomplete zeta function and 

y p (u,y):=-log(l-uj) + Y^?^(-iy(l-aj)- r - f C 2p+1 (x) x - 2 - 2 P dx. 

r=l T Jl ~ u 

Proposition 4.4. Let lu — 0,1/2. Then 



*p(w,«) = Hm(C p (w,i/;s) - !/(*-!)), 



C P (w,y;-n) =-— =C(-n), n = 0, l,...,2p, 
F n + 1 

C(s) = lim C p (w,y;s), Res + 2p>0, 



* p (w,y)-7= / C 2p +i{x)x 2 2p dx, 
Jy 

7 = lim %(<j,y). 



Proof. The second relation follows from g3 Eq. 2.8(13)], S 2fe +i(w) = for uj = 0, 1/2 and fc > 1 
and [TJ Eq. 23.2.15]. The representations and therefore the limit relations for £(s) and 7 follow 
from Proposition ^. II □ 



Proof of Theorem \2.3l Let / be admissible in the sense of Definition 12.21 In the representation 
(|4.4p we can write the integral as follows: 



o Hri/2 2 /" |r|/2 2 /" |r|/2 

m/. /w *-m/. s * )dI+ m/ </-«.><*>* 

2 * /•l r l/ 2 2 f |r|/2 2 

= ]f[ 2-] _/ X_;5 " + ]f[ i (f-S q ){x)dx~ — {f~S q )(x)dx 



V f y 



I I n=Q s « I I JO 

Defining 

rj /"^i — U>;|1 |/JV 

$n p (/-^;7V):=-— iV 2 y (/-5 g )(x)dx-^(r,/-5 g ;iV)+^ p (r,/-5 g ;7V), 



-(1— u>)|r|/iv 



formula (14.41) becomes (in condensed notation) 
2 9 a 1 " 8 ™ 

M(f; N) = V f N 2 - —N 2 £ ^y^- - ^,(5,; JV) + B p (f; N) + K p (S q - N) + D\ p (f - S q ; N) 



|r| 

1 1 n=0 

= Vf n 2 + J2 ^ I y\ n2 T~i ~ • A *( x ~' n ; N) + ; ^)} + W-> N ) + - S i : 



n=0 
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Furthermore, using (|4.5a|) , (14.5c|) and Definition 14.31 we can write the expression in curly brackets 
above as follows: 



2 , a 1 -*- 
— N — 
iri S r , 



2p 



- - A p (x~ s - ; N) + K p (x~ s " ; N) = Atf3±Jl _ ^2 £ 



B r (u) 



\T\/N) r {t-^} {r 



2N 



(\T\/N) 



2p+l f N/2 



(2p + l)l h 



C 2p+l (x) {t— } 



(2p+l) 



dx 



= — N 2 (\T\/N) 



t=x\r\/N 



( S ")2p+1 



N/2 



Hence, we arrive at the formula 



C^+t^x-^-^dx^ —AT 2 (|r] /iV) 1 " 8 " C P (^, iV/2; 



M(f;N) = V f N 2 + J2 



" 2C p (u,,N/2; Sn ) 



N 1+s " +B p (f;N)+m p (f-S q ;N). 



n=0 



For m p (TJ;N) defined by (JHHJ) we have 



J-^ 2(JK/2,N/2:s n )-2((sn) ^ 
(4.7) V\ p (T, /; N) = J2 a n ^ ' |S „ ^iV 1+s " + m p (f - S q ; N). 

n=0 ' ' 

Furthermore, it follows from Lemma g21 that «H P (/ - S q ;N) = OiN 1 -^ 3 ") + OiN 1 ^) if 2p 
6 - Res, and <K p (/ - S q ;N) = 0(N 1 ~ 5+S ") + ©(iV 1 -^ log AT) if 2p = 6 - Rcs q . Finally, using 
(|4.6[) and Proposition 14.41 we obtain the estimate 



* C p W2,iV/2, S »)-C(^) Arl+Sr 

71=0 



, r ,-„ 



< 2(iV/2) 1 - 2p ^ 



71=0 



_B2p_ v 2p + s n 
an (2p)\ [Snh P2p + Re Sn 



(|r|/2)" 



- Rc St, 



Note that, whenever s„ = — k for some k = 0,1,... ,2p, then the corresponding terms on both 
sides of the estimate above are not present. Also, from Definition ^. 2l it follows that 2p + Res„ > 
for n = 0, . . . ,q — 1 and that either Re s q + 2p > or s q — —2p. In either case the sum on 
the left-hand side above is of order 0(N 1 ^ 2p ). Hence, we have from gj]) that m p (TJ;N) = 
0(N 1 ~ s+s ->) +0(N 1 ~ 2 p) if 2p ^ 6 - Res q and 9t p (r, /; N) = OiN 1 -^ 3 ") + 0(N 1 ~ 2p log A^) if 
2p = S -Res q . □ 



Proof of Theorem \2.J\ Proceeding as in the proof of Theorem 12.31 the remainder term now takes 
the form 



*p(r, /; N ) = w\ N2a «' N l 2 ) - 7) + E a ™ 



2C p («/2,JV/2, Sn )-2C(an) Arl+a „ 



n=0, 



— AT — 



(i-a>)|r|/jv 



(/ - S q )(y) dy - A P (T, f - S q ; N) + TZ p (T, f - S q ;N). 



Using Lemma |4~21 Proposition 14.41 and the inequality 



Ajv%, (-1^/2^/2) - 7 ) 



< 4— \a q ,B 2p \ (N/2) 



l-2p 



we get the estimate 9t p (I\ /; AT) = 0(N 1 ~ 2p ) + 0(N l - s+s «) if 2p f S - Rc s q and m p {T, /; AT) = 
OiN 1 - 2 ^ logN) if 2p = 6 -Res q . ' □ 

Next, we prove the results related to particular types of kernel functions. 
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Proof of Theorem \2.5\ The Laplace transform f(x):= e~ xt d/i(t) of a signed measure fi on 
[0, oo ) satisfying t m d\/i\(t) < oo for every m = 0, 1, 2, . . . has derivatives of all orders on (0, oo). 
For q a positive integer let S q (x) be defined by S q (x):= J2n=o ^(~ x ) n - ^ 0T ever y < ra < g we 
can write 

Jo ^ m (n-m)\ 

where, using a finite section of the Taylor series expansion of h(x) = e~ xt with integral remainder 
term, we have that 



poo I 1 m ( j.\n 

(/ - S q )< m \x) = f^\x) - S^(x) = (-l) m / { e' xt - V 



n=0 



t m dn(t) 



_ (~l) g+1 
(q-m)\ y 

For x > we have the following bound 



j jj^ e- ut (x - u) q - m du\ t q+1 dfi{t), x>0. 



T g+l — m 

(f-S q )( m \x) < _ J m = 0,l,. 



,q. 



Since S q (x) = for all x, it is immediate that the last estimate also holds for m = q + 1. It 
follows that / is admissible in the sense of Definition 12.21 with q = 2p, 6 = 1. The result follows 
from Theorem 12. 31 after observing that 

o f|r|/2 2 f |r|/2 f°° 

V f (r) = ^ /(*) <fc = p / n J e- xt dix{t) dx. 

□ 



In the case that / is a completely monotonic function on (0, oo) (that is, /i is a positive measure), 
it is possible to improve the estimate for TZ p (T, /; N) in (|4.5c[) . 

Proof of Theorem \2.6\ Let / be analytic in a disc with radius |T|/2 + e (e > 0) centered at the 
origin. Then f{z) = X^^Lo a n z ™ f° r \ z \ < |r|/2 + e and / is admissible in the sense of Definition ^. 21 
for any positive integers p and q — 2p, where S2 P (z) — X)n=o a « z " an< ^ 5 = 1. The asymptotic 
expansion follows from Theorem 12.31 on observing that with s n — —n (n = 0, . . . , 2p) 1 one has 
2 2p |r| /2 pi/2 2 .|r|/2 

W = pE a «{ xndx + W\Jo U- s ^& dx =ir\J o fW dx - 

Moreover, since s q = -2p and 8 = 1, it follows that «H p (r, /; iV) = e> p ,|r|. / (iV 1 ~ 2p ) as iV -> oo. □ 

Proof of Theorem \2.8l Suppose / has a pole of integer order K > 1 at zero and is analytic in 
the annulus < \z\ < \T\/2 + e (e > 0) with series expansion f(z) — 53nL— k a « z ™- Then / is 
admissible in the sense of Definition 12.21 for any positive integers p and q — 2p with S2 P (z) — 
Yl 2 n=-K anZn an d ^ = 1. In the case (i) Theorem 12.31 is applied and in the case (ii) Theorem 12.41 
is applied. The expressions for Vf(T) follow from termwise integration in (|2.4[) and (|2.5p . Since 
l-S + s q = -2p, the remainder terms are 9\ P (T, /; N) = O p ^ T \ i f(N 1 ~ 2p ) as N -> oo. □ 

Proof of Examples \2.9\ and \2.l0[ If / has an essential singularity at and is analytic in the annulus 
< \z\ < \T\/2 + e (e > 0), then for positive integers p 

2p oo 

f(z) = S 2p (z) + F 2p (z), S 2p (z):= a n z n , F 2p (z):= £ a n z n = 0(z 2p+1 ) as z -> 0. 

n— — oo n— 2p+l 

Clearly, the function /(z) satisfies item (i) of Definition 12.21 and both functions f(z) and S 2p (z) 
satisfy an extended version of item (ii) of Definition 12.21 suitable for an infinite series S 2p (z). Since 
termwise integration and differentiation of S 2p (z) are justified by the theory for Laurent series, 
Theorems 12.31 and 12.41 can be extended for such kernel functions /. In this case all formulas in 
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Theorems 12.31 and 12.41 still hold provided the index n starts with — oo. In particular, we note that 
the infinite series ^n=-oo n#-i a ™ C(~ n ) 1^1" N 1-n appearing in the asymptotics of M(T, /; N) 
converges for every TV, since C( m ) < C(2) for all integers m > 2. 
Example 12.91 follows from the extended version of Theorem 12.41 

To justify Example 12.101 let A be a zero of the Bessel function J_ x . The extended version of 
Theorem 12.31 with a n — J n (A) gives that for positive integers p > 2 and m > 2 

oo 

M(T,f;N) = V f (T)N 2 + 2 ]T J-n(A) C( n ) |rp" N 1+n + B P (T, f; N) + 0{N 1 ^ 2p ) 



n——2p, 



2N £ J_ n (A) C(n)(N/ \T\r + 2 ]T J_„(A) C(n) \T\~ n iV 1+ " + 7/ (T) iV 2 + |r| B 2 £)f>0) 



n—m 7i—2 
2p 



2" v 2 



+ 2£ J fe (A) C(-fe) |F| fe iV 1 ^ + £ ^2fc^-/0»-D(|r| /2)iV 2 - 2 » + O(N^), 

k=2 n=1 ( 2n )- l r l 

where 

V (v) 2 V t r.v / |r|/2)1+ " 

1 1 n— — oo. 

In the above we used the relation (12. 2j) . Observe that £(— fc) = for fc = 2, 4, 6, □ 



Proof of Theorem \2.11\ The asymptotics and the remainder estimates follow from Theorem [2? 
on observing that ff{x) has derivatives of all orders in (0, |T|/2 + e), S q (x) = 2n=o a n% n ~ s , and 
5=1. The constraints on s q — s — q imply that the positive integers q,p and s € C satisfy 
q — 2p < Res < 2 + g or s = q — 2p. For < s < 1 we have (see (|2.4jl ) 

2 r' r|/2 ^ 2 ~ (|r|/2) 1+ — 



9 2 fir 



and the right-hand side as a function of s is analytic in C except for poles at s = 1 + n (n = 
0,1,2,...) provided a„ ^ 0. □ 

Using the same method of proof as in 20J for the Hurwitz zeta function, we obtain the following 
two propositions, which will be used in the proofs of Theorems 13.21 and I 



Proposition 4.5. Let q > 1 and a — 1/2 or a = 1. For a; > and s E C wit/i s ^ 1 and 
Re s + 2q + 1 > the Hurwitz zeta function defined as £(s, a):= J2kLo(k + a )~ s / or Rc s > 1 and 
o^O, — 1, —2, . . . has the following representation 

l — s Q f \ 

C(s, z + a) = f— j - (a) aT s + ^ 2 " " («) 2 n-i a;1 ~ 8 ~ 2 " + P?( s > 

n— 1 ^ '' 

The remainder term is given by 

Pq {s, x, a) = -L P +l °° EHgjrVMO c(s + rfw = a (ar i-Re -2 9) 

2tt« J 7 „-,:oo T(s) 
os iV — > oo, where —1 — Re s — 2q < ^ q < — Re s — 2q. 

By the well-known relation log[F(a; + a)/\/27r] = C( s j x + a )|s=o one obtains the next result 
from Proposition ^. 51 



Proposition 4.6. Let q > 1 emd a = 1/2 or a = 1. For x > 

log = (x + a - 1/2) logx-x + Y, <2n-l)2n xl ~ 2n + P ^ ^ 

V 71— 1 ^ ' 
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The remainder term is given by 



Pi( x > a ) = I T(w) C(w, a)x w dw = Ogix- 1 - 2 ") 



2ni 

as N — > oo, where — 1 — 2q < j q < —2q. 



In the proofs of Theorems 13.21 and 13.31 we make use of the observation that for N — 2M + K 
with M > 1 and k — 0, 1 formula (|2.ip simplifies to 

2 LA72J 1 1 - 

(4-8) M S (T;N) = — -N 1+s V — - „ , , 3 JV. 

I r l nrf fcs (l r l/2) 

which involves the generalized harmonic numbers i?„ -=X]fe=i ^ — 



Proof of Theorem \3.SX Differentiating (|4.8[) with respect to s and taking the limit s — > yields 

AT 

X log (r ; Af) = A" (N — k) log — - 27Vlo g r(LiV/2j + 1) - (1 - K)N\og(N/2). 

The asymptotic expansion of the theorem now follows by applying Proposition 14.61 with x — N/2, 
a = (2 - k)/2. Note that B 2n (a) = B 2n (l - k/2) = B 2u {k/2). □ 



Proof of Theorem \3.3[ Starting with Theorem 12 - 31 we obtain an asymptotic formula of the form 
(pH]) but with error estimate 0( AT 1-2 9)0 On the other hand, substitution of the identity 
Y^k—i k s = C( s ) — C( s , n + 1) m to l|4.8p gives the exact formula 

M S (T:N) = - -^ 1+S C(*, L^/2J + 1) - 

Then the asymptotic relation (|3.2p with error term of order 0(N~ 2q ) follows by applying Propo- 
sition [43] with x = N/2, a = (2 — «)/2. This expansion holds for s with Res + 2q + 1 > 0, 
g>l. □ 

Proof of Proposition \3.4] Using Jacob Bernoulli's famous closed form summation formula ([TJ 
Eq. (23.1.4)]) V + 2P + ■ ■ ■ + nP = (n + 1) - B p+1 )/(p + 1) in (TCT) one gets 

M-JT; N) = 2 Tf ^P+i(( iV + K )/ 2 )-- B P+i jV i- P + (!_*) (|r| /2f N. 

p+l 

Application of the addition theorem for Bernoulli polynomials (see [TJ Eq. (23.1.7)]) yields the 
result. □ 

Proof of Theorem \3.5[ An asymptotic formula with error estimate 0(N 1 ~ 2p ) follows from Theo- 
rem [2j4j see also the second remark after Theorem 12.41 However, by substituting into (|4.8p with 
lu = k/2 (k = 0, 1) the following asymptotic expansions 

/a q\ „ v^l , f , A , fi iM -B2fcM/(2fc) i? 2g+2 (a;)/(2 g + 2) 
(4.9) H n = ^- = log(n + w) + 7 - — - - 2^ ; , ^ ± , ~ ^+2 » 



where < q ^,K < 1 and collecting terms we get the asymptotic formula (|3.4p with improved 
error estimate. The plus sign in (14. 9p is taken if w = 1/2 and the negative sign corresponds to 
oj = 0. We remark that the representation (|4.9p is given in [11] if a; = 1/2 and can be obtained 
as an application of the Euler-MacLaurin summation formula if co = (see, for example, [3]). We 
leave the details to the reader. □ 



ttlf Res = 



— 2g and s ^ 2q, then a factor logN must be included. 
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